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1. Introduction. It is well-known that a number of steady state prob- 
lems in fluid mechanics involving systems of nonlinear partial differential equa- 
tions can be reduced to the problem of solving a single operator equation of the 
form 

( 1 . 1 ) v + XAv -I- XB(v) = 0 , v E H , X G R 1 , 

where H is an appropriate (real or complex) Hilbert space. Here X is a typi- 
cal "load’' parameter, e.g., the Reynolds number, A is a linear operator and B 
is a quadratic operator generated by a bilinear form. In this setting many 
bifurcation and stability results for problems in fluid mechanics have been 
obtained and the reader is referred to [ 1 ; 10 ; 11 ; 21 ] and the bibliographies 
therein for a detailed account of such results. 

In fact, there may be considerably more structure in a nonlinear stability 
problem in fluid mechanics than that implied by an operator equation such as 

(1.1) . As shown in a recent series of papers by the authors [12;13;19], a "struc- 
ture" parameter, say 7 , also may be present so that equation ( 1 . 1 ) may be actu- 
ally of the form 

(*) u - X(L-7M)u - F(u) - tG(u) = 0, u EH, X 6 R 1 , 7 G R 1 , 

where H is again an appropriate Hilbert space, L and M are linear opera- 
tors, and F and G are generated by bilinear operators. Equations of the 
form (*) are derived in [12;13] for Benard-type convection problems and in [19] 
for the Taylor problem. It is the purpose of the present paper to describe a 
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setting in which it is possible to determine a complete set of bifurcation equa- 
tions for an operator equation of the form (*) by using the structure parameter 
7 as an "amplitude" parameter rather than regarding 7 as merely a constant 
to be incorporated into the operators M and G in (*). Complete bifurcation 
diagrams are obtained by such an approach because, e.g., the "reduced" bifur- 
cation equations contain linear as well as both quadratic and cubic terms. This 
is in contrast to many other approaches in bifurcation theory based upon split- 
ting methods (e.g., see [ 10 ; 11 ; 18; 21 ]) in which either linear and quadratic or 
linear and cubic terms appear in the reduced equations. In this sense the spirit 
of the approach presented here is somewhat like that of singularity theory as in 
Golubitsky & Schaeffer [7], however, our approach also has points in common 
with those in Busse [ 2 ] and Chow, Hale & Mallet-Paret [4; 5]. The present 
paper emphasizes the identification and utilization of structure parameters in 
fluid mechanics but the approach presented here applies equally well to general 
classes of bifurcation problems provided that such problems have appropriate 
"higher order" terms (e.g., see the general class of variational problems studied 
in [14]). 

The following hypotheses are chosen to illustrate our approach with a 
minimum number of technical difficulties. We consider equation (*) under the 
following hypotheses in which H is a real Hilbert space with norm II II gen- 
erated by an inner product ( , ). 
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(Hi) The linear operator L: H — ► II is selfadjoint and compact, and, if 
Ho clenotes the smallest positive characteristic value of L, then the dimension 
of the null space N = A r (I— /u 0 L) of L is n = 1 and N is spanned by ip 0 
with lh/> 0 ll = 1 . 

(H 2 ) The linear operator M: II — ► II is compact and satisfies 

M: N — ► n\ where denotes the orthogonal complement of N in H. 

(H3) The critical characteristic value, \ c = \{l), of L 7 s L — 7 M, i.e., 
the positive characteristic value of L 7 of least magnitude, is of the form 

(1.2) X c (7) = Ho ~ HoW + A(t), hi < Tu 

where A is real and analytic and satisfies A( 7 ) = 0(t 3 ) as 7 — +■ 0. Here 
b = (MKMVb, 0 o) where Iv: denotes the inverse of the restriction of 

(I— PqL) to N . 

fH4l The nonlinear operators F(w) == <J»(w,w) and G(w) = T(w,w) are 
generated by bounded bilinear operators H X H — *■ H and 

T: H X H —*■ H, where, in addition, N X N —+ N . 

The assumption in (Hi) that dim N = 1 is made merely for the sake of 
convenience; more general settings with dimA r > 1 may be treated along the 
lines in [12; 13; 10], however, the actual solution of the bifurcation equations 
may be, of course, much more difficult in the case dimA r > 1 . The assumption 
(1.2) in (H3) is the "natural" expansion for X c and also is made for the sake of 
convenience. It, or a suitable modification, holds in many linear problems in 
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fluid mechanics even if M is not symmetric and dimiV > 1 (e.g., see [3; 8; 12, 
13; 19]); if M is symmetric, then (H3) may be omitted. The modifications 
required when H is a complex Hilbert space are somewhat more involved but 
may be formulated as in [12, §4] in terms of real operators and group represen- 
tations (see also [20; 21]). Finally, if one wishes to carry out a linear stability 
analysis of solutions of (*) in the case dimiV > 1, one may proceed as in [12, 
§5] (see also [16]). 

Remark 1.1. In many problems in fluid mechanics the bilinear operator 
$ in (H4) will be generated by the nonlinear term (u-V)u in the Navier-Stokes 
equations. If so, then $ satisfies the additional condition 

(1.3) ($(u,v),w) = -(<f>(u,w),v), u,v,w € H . 

The condition (1.3) plays a role in determining the actual form of the bifurca- 
tion equations associated with (*) but it is not required in the derivation of the 
bifurcation equations. Thus, condition (1.3) is not included as part of (H4), 
however, it is used in the application to rotating Couette-Poiseuille flow in Sec- 
tion 4 and also in the applications to Benard-type convection problems in [12; 
13] and to the Taylor problem in [19]. 

The outline of the paper is as follows. In Section 2 we make use of 7 as 
an amplitude parameter and derive the bifurcation equation associated with (*) 
under hypotheses (Hi) through (H4). In Section 3 the bifurcation equation is 
solved under an additional invariance assumption on the remainder term and, 
for each fixed 7 sufficiently small, a complete bifurcation diagram is obtained. 
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In Section 4 we use the above approach to study rotating Couette-Poiseuille 
channel flow. We show, in particular, that the superposition of a Poiseuille flow 
on a rotating Couette channel flow is, in general, destabilizing. This type of 
result was conjectured in [17] for non-rotating combined Couette-Poiseuille flow 
on the basis of numerical calculations for the linearized Navier-Stokes equa- 


tions. 
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2. The bifurcation equation. The bifurcation equation associated with 
(*) is derived here using standard splitting methods except that the structure 
parameter 7 plays the role of an amplitude parameter. 

Let P: H —*■ N denote the orthogonal projection of H onto N and 
let S = I — P: H —*■ N denote the orthogonal projection of H onto N. For 
7 sufficiently small, we shall seek solutions of (*) of the form 

(2.1) u = 7 (i/H- 7 ^), X - Mo - /vrVo b ~ T )i 

where ip € N, ^ G N^" and r G K 1 are to be determined. Note that if 
t = t 0 + 0(7) and if X c = X c (7) is defined as in (H 3 ), then 
X = X c + /UoTqT 2 + 0 ( 7 *). Thus, for 7 sufficiently small, a solution of (*) of the 
form (2.1) is subcritical if r 0 < 0 and supercritical if r 0 > 0. 

Substituting (2.1) into equation (*), using the projection P onto and 
S onto N, and making use of (H2) and (H 4 ), one obtains the following equa- 
tions in N and N: 

( 2 . 2 a) 0 = (I— /u 0 L)'I' + n 0 M ip — F (ip) 

+ iP[/J, 0 My—$(ip,ty)—$(ty,ip)—G(ip)] 

+ l s PMr-|Uo b ) M '/' + b-r)L# 

- f(») - r (**) - r(^»] 

+ T 3 PMr-,u 0 2 b)M'l' - GOJ-)], 

(2.2b) 0 = (ftj 2 b-r)V’ + 

- TSlF(*)+r(*«)+r(W)| + T 2 SMr- rt ?b)M'J'-G(1<)]. 
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Since K = [(I— /^L)!^] 1 is bounded on iV"7 given p > 0 there exists 7 0 > 0 
such that, if (ip, r ) E N X R 1 with |r| < p and WipW < p, then by the implicit 
function theorem (2.2a) can be solved uniquely for = 'l>(ip,T,'i) in N pro- 
vided that 1^1 < 7 0 . In fact, V& is analytic and of the form 


(2.3) * = —p 0 KMip + KF(</>) + 7*1, 

where ^ ^ l (ip,T,'i) E N 1 is bounded with bound depending only on p. It is 

important here that p can be arbitrarily large provided that 7 0 is chosen 
sufficiently small. Subbing ^ into (2.2b), taking the inner product with ip 0 , 
and making use of the definition of b in (H3), one obtains, for ip of the form 
ip = (3ip 0 with fS E R 1 , the bifurcation equation associated with (*): 

(2.4) 0 = -t( 3 + (P[p 0 (MKF(ipQ),ip 0 ) -f p Q (<t>(ip 0 ,KMip 0 ),ip 0 ) 

+ p 0 ($(KMip 0 ,ip 0 ),ip 0 ) - (G(4),-0 O )| 

- l^[(^(ip 0 ,KF(ip 0 )),ip 0 ) + ($(KF(ip 0 ),ip 0 ),ip 0 )} + r (P,T,i) 

= —t/3 + a/3 2 + CjS 3 + r(P,T,-y) = Z>(/?,r,7). 

Here, for |r| < p, \/3\ < p and I7I < 7 0 , the remainder term r is given by 

(2.5) r(/?,r,7 ) = (7S{m 0 M^ 1 -4>(i/>,^ 1 ) - ^^ip) 

- F(tf) - T(iP,V) - T(%iP) 

+ 7K( r “^o b ) M ^ “ 


where ip = plip 0 and ^ = '^ l (f3ip 0 ,r,'f). Moreover, r is analytic in (/3,t, 7) 
and, for some r 0 depending only on p, satisfies 
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(2.6) |r(/?,r,7)l < r oM> 1^1 < P, M < P< hi < 70- 

Note that if, for fixed 7 satisfying hi < 7o> (/?*,f*,7) * 1S a solution of 
(2.4) with \r*\<p and \/3*\<p, then (v*, X*) G H X R 1 given by 

(2.7) v* = 7(Wo + 7'J'(/?*Vb> r *i7), X* = ^0 ~ Mo7 2 (^o b - 7 *) 
is a solution of equation (*). 

In the next section we solve the bifurcation equation (2.4) under an addi- 


tional invariance assumption on the remainder term r. 
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3. Solutions of equation (*). In many problems of the type considered 
here the remainder term r in (2.5) may satisfy additional properties. E.g., 
because of invariance and symmetry considerations, r may have the form 

(3.1) r{P,T,i) = /&(/?, r, 7 ), 

where s is analytic and s and its partials with respect to ft and r are uni- 
formly of order 0 ( 7 ) as 7 — *■ 0 ; one has, e.g., 

(3.2) kAvy)l < s 0 hi, \P\ < P, M < P, hi < 7 o, 

for some s 0 depending only on p. In such problems the bifurcation equation 
(2.4) may be factored and replaced by 

(3.3) 0 = -T + a/? + cfi 2 + s(P,T , 7 ) = E{P,t, 7 ). 

It is this type of "factoring” that is also the key to the solution of more difficult 
problems with dim N > 1 (e.g., see [13; 15; 19]). Since E is a mapping of a 
neighborhood of (0,0,0) £ 1K° into JK% it is natural to seeK solutions of ( 0 . 0 ) 
for 7 near 7 = 0, hence solutions (v*, X*) of (*) in the form (2.7), by means of 
the implicit function theorem. We have, e.g., the following result. 

Theorem 3.1. Given p > 0 there exists 7 0 > 0 such that for hi < 7o 
equation (3.3) has a solution 

(3.4) r = t(/?, 7 ) = a/? + eft 2 + r^/3, 7 ) 


that is bounded, analytic and unique in 
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(3.5) B = |r— a/?— c^ 2 1 < kpy|, \0\ < p, HI < 7o}, 

where the constant k depends only on p, and Tj is 0 ( 7 ) as 7 — ► 0 , uni- 
formly for |/?| < p. For each fixed 7 satisfying [ 7 ! < 7 0 , the corresponding 
nontrivial solution branch (v*(/?), X*(/?)) of (*), \/3\ < p, has the form (2.7) with 
f3* — 0 and 7* given by (3.4). 

Proof. Let /? 0 satisfying \(3 Q \ < p be given and set r 0 = a /? 0 -f c/3q. It 

d .E 

follows from (3.2) that {/3q,t 0 ,0) is a solution of (3.3) at which = —1. 

Thus, by the implicit function theorem, (3.3) has a solution r == t(/ 3,7) that 
is bounded, analytic and unique in a neighborhood of (/?,r, 7 ) = (/? 0 ,r 0 , 0 ) with 
r{@ 0 ,0) = r 0 . Since 

(3.6) |r-a/?-c^| = |s(/?,r, 7 )| < s 0 HI < s 0 7o> 

a finite number of such neighborhoods cover B, provided that 7 0 is sufficiently 
small. For fixed 7 satisfying I 7 I < 7 0 , the form of the corresponding solution 
branch (v*(/?), X*(/?)) £. H X R 1 is an immediate consequence of the above con- 
struction. 

Remark 3.1. One can show also that, for each fixed 7 satisfying 
HI < 7 oj if c > 0, then the branch of solutions, Q, of (3.3) determined by (3.4) 
has a unique turning point in B at /3 = /? T ( 7 ), where /? T is of the form 

(3.7) Pr(l) = — (a/2c) + ^( 7 ) 

with P \(l) = 0 ( 7 ) as 7 — *> 0 . The fact that (/?, 7 ) =0 in (3.4) at a 
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unique point in B follows from another application of the implicit function 
theorem; one considers the equation 


(3.8) 0 = 4L = a + 2c/* + 4L r,(, 3,i) 


c)t 

and uses that, at '7 = 0, — — = 0 only at (3 = —a /2c. 

op 


It can be shown (e.g., see [12, §5; 16; 21] that, for fixed 7, the linearized 

stability at points (v*, X*) along the solution branch of (*) in Theorem 3.1 is 

dD 

determined by the sign of — — along the branch Q. Since, along Q, (2.4) 

dp 


holds and 


(3.9) |2- - /?(a+2c/?+||), 

one sees that, for 7 0 sufficiently small, the stability of ( v*(P ), \*(/?)) is 
indeterminate only at p = 0 and at the unique turning point on Q given by 
(3.7). Thus, we have the following corollary to Theorem 3.1. For the statement 
of the corollary we may assume that a < 0; if a > 0, one replaces P by — P 
in the given intervals. 


Corollary 3.1. If a < 0 and c > 0 then, for 7 0 sufficiently small, the 
solution branch (v*(/3), X*(/3)) of (*) obtained in Theorem 3.1 is stable for 
— p < P < 0 and /?x(t) < P < P and unstable for 0 < P < /3 T ( 7). 

We wish to emphasize that, for each fixed 7 sufficiently small, the solu- 
tion branches of the bifurcation equation (2.4) obtained in the above discussion 
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are "global" in the {/3,t) plane, however, the corresponding solution branch 
(v*, X*) of (*) may, of course, be "small" because of the amplitude factor -y 
in (2.7). 

The above results are formulated only to illustrate the approach when 
dimiV = 1. Complete solutions of the bifurcation equations and the 
corresponding stability results may be, of course, much more difficult to obtain 
in problems where dim N > 1 because in such problems, e.g., there may be 
points of secondary bifurcation at which the null space of the appropriate 
linearized operator has dimension greater than one, (see the discussion of hexag- 
onal solutions in [13]). 
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4. Rotating Couette-Poiseuille channel flow. In this section we study 
rotating Couette-Poiseuille channel flow and show that, in general, the superpo- 
sition of a Poiseuille flow on a rotating Couette channel flow is destabilizing. In 
the case of non-rotating Couette channel flow this result was conjectured in [17] 
on the basis of numerical calculations for the linearized Navier-Stokes equa- 
tions. In the nonlinear analysis presented here, it is crucial that there are 
Coriolis effects present in the problem so that the swirl-like parameter S 
defined in (4.2) is positive. 


We consider viscous incompressible flow in a rotating infinite-channel of 
width l. The non-dimensional Navier-Stokes equations in a rectangular coordi- 
nate system rotating about the z-axis with constant angular speed, fi, are given 
by (e.g., see [9, p. 163]) 


(4.10).| R e *Av — Vp' — (v • V)v + 2S 


Vo 

— V, 


= 0 , 


(x,y,z) G C, 


^V-y = 0 


where 


C = {(x,y,z): — oo < x < oo, 


< y < -oo < z < oo}, 


and 


v = (v 1 ,v 2 ,v 3 ). Here the Coriolis acceleration terms are determined by the 
square brackets, the rectangular coordinates are scaled by l, the velocity com- 
ponents are scaled by U c , where U c denotes the maximum velocity of pure 
Couette flow at y = l/2 in the direction of the x-axis with no Poiseuille flow 
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present (e.g., see [9, p. 179]), S is the swirl-like parameter 

(4.2) s = m/u c , 

and R e is the Reynolds number 

(4.3) R e = lU c /v, 

where u is the kinematic viscosity. We assume throughout that S is fixed and 

0 <S < -7. 

2 

For the basic unperturbed flow we take a pressure, P, and a combination 
Couette-Poiseuille velocity distribution U = (U x (y), 0, 0), where 

(4.4) U,(y) = (y+i) + | (1— ly 2 ), - < y < i , 

(4.5) 6 - U p /U c . 

Here U p denotes the maximum velocity of pure Poiseuille flow in the direc- 
tion of the x-axis with no Couette flow present (e.g., see [9, p. 66]). 

We shall seek solutions of (4.1) that are perturbations of U and P, and 
are independent of x. Setting v = U + w, p' = P -f q in (4.1) and assuming 
w = w(y,z) and q = q(y,z), one obtains the disturbance equations 

f W 2 

R e ^w — Vq + 2S — w } — w-VU — w-Vw = 0, 

0 

Y'w - 0, — j < y < j, -oo < z < oo, 

d d 

where V = (0, The boundary conditions here are w = 0 on 
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y = ±-^-. To introduce the appropriate parameters we set 
£ 


(4.6) X = [2S(l-2S)] 1/2 R e , 0 < S < 

2 


(4.7) 


Wj = X 1 (2S— l)u 1? w 2 = X 1 [2S(1 — 2S)] 1 / 2 u 2 

w 3 = X- 1 [2S(l-2S)] 1 / 2 u 3 , q = X- 1 R e " 1 [2S(l-2S)] 1 / 2 p, 


and define the structure parameter, -y, as 


(4.8) 1 = 


— = Hp o < s < 1 

1-2S U c (l— 2S) ’ ^ 2 


The disturbance equations for u = (u 1? u 2 , u 3 ) and p become 


(4.9a) 

(4.9b) 

(4.9c) 

(4.9d) 


Auj + X(l— 2-yy)u 2 - 


5u t <9u! 


U2 ^r + ' 13 az 


2 9y 


+ Xuj — 

( Pn 


.9 <9uf> 

U 2 -ZT + u 3 


du 2 

dy 


Au3 “af- 


A. 


U2 1T + U3 ^ 


dz 


= o 


= 0 


= 0 


<9u 2 <9u 3 i i 

■^r + "5T“ 0 ’ - l <y< l’ -° 0<z<0 °- 


with the boundary conditions 


(4.10) u i = u 2 = u 3 = 0 on y = ±-^ . 

The equations in (4.9) are^ closely related to those for the generalized Benard 
problem studied in [12; 13]. In carrying out the analysis below for rotating 
Couette-Poiseuille flow we shall make repeated use of this relationship. 
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We next introduce an appropriate Hilbert space setting in which to seek 
solutions of (4.9) that are periodic in z. Given a positive number a (to be 
specified below in (4.22)) we set 

(4.11) R = {(y,z): < y < ~ and 0 < z < — }. 

2 2 a 

The Hilbert space, H, used throughout this section is defined as the closure of 
the set {v = ( v i*v 2 ,v 3 ): v is smooth, periodic in z with period and van- 

Oi 

ishing in a neighborhood of |y| = — with V-v = 0} in the norm, II II, associ- 

Zi 

ated with the inner product 

(4H2) (v,w) = / £ Yvj’Ywj • 

R j=i 

Here and in the sequel, whenever possible, the vector notation v is suppressed 
when dealing with elements of H. 

To formulate the problem as an operator equation in H, we take the 
scalar product of (4.9a) through (4.9c) with w£if, and use (4.10) and integra- 
tion by parts to obtain 

(4.13) (u,w) — X(L 7 u,w) — (F(u),w) = 0. 

Here, for each 7 E R 1 , the linear operator L 7 : H — ► H and the quadratic 
operator F: H —*■ H are given by 

(4.14) L 7 = L — 7M, 

(4.15) F(u) = $(u,u), u EH, 
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where the operators L: H — ► H and M: H — ► H are defined (weakly) by 

(4.16) (Lv ,w ) = / (v 2 w x + v 1 w 2 ), 

n 

(4.17) (Mv,w) = 2 f yv 2 W!, v,w, £ H, 

n 

and the bilinear operator H X H — ► H is defined (weakly) by 

(4.18) (4>(u,v),w) = -/ [(u-V)v]-w = -/ (u 2 ~ + u 3 — -)-w, 

Since w in (4.13) is an arbitrary element of H, one obtains an operator equa- 
tion of the form (*) in Section 1, namely 

(f) u— \L 7 u — F(u) =0, u £ H, X £ Ft 1 , 7 £ Ft 1 . 

Standard regularity methods (e.g., see [10; 11]) can now be used to show that 
the problems of finding classical solutions of (4.9), (4.10) and solutions of ( f ) in 
H are equivalent. 

We require the following facts about the linear problem associated witn t TJ 
when 7 = 0, namely the problem 

(4.19) u — /uLu = 0, u £ H, fj, £ Ft 1 . 

The linear problem (4.19) is equivalent to the classical problem, for smooth u 
and p periodic with period 27 r/a in z, obtained by setting 7 = 0 and omit- 
ting the nonlinear terms in (4.9). It is sufficient to consider the solutions of 

(4.19) given by (see also [12, (2.1 l)ff] and [10, §3]) 

(4.20a) u = (<^ 1 (y)cos<72, 0 2 (y)cosoz, ^ 3 (y)sinOT), 
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(4.20b) p = —a 1 cosazD 2 (f) 3 , 

(4.20c) </> 3 = -cT l <f> 2 , 

O d 2 d 

■where D 2 = — - — a 2 , a prime denotes , and and (f) 2 satisfy 

dy J dy 

(4.21a) D V 2 - fid 2 <j>i = 0, 

(4.21b) D 2 <j>i + fi(j) 2 = 0, 

(4.21c) <£i = <fh = (f>2 = 0 at y = ±~ • 

(The solutions of (4.19) with cosot and sinOT interchanged lead to flows that 
differ only in orientation.) One can show for cr > 0 (e.g., see [8]) that the 
eigenvalue problem (4.21) has a countable number of positive, simple eigen- 
values, 0 < fiy(o) < fi 2 (o) < • • * , depending continuously on cr. Moreover, 
— * oo as either cr — ► 0 + or cr — ► oo. Consequently, fi\{cr) assumes an 
absolute minimum at some <x 0 > 0. We assume that <r 0 is unique so that 
t-h{°) > M i(°b) */ ^ <r 0 . (This property is suggested by numerical calculations 

(e.g., see [3, §15(b)] and [6, §10]) and is usually assumed in such problems.) For 
some given integer p 0 > 1 we now choose a so that 

(4.22) <t q = p 0 o; 

and use this a to define the basic Hilbert space H of this section. 

We may now determine a complete solution of the linear problem (4.19). 
Since cosot and sinOT in (4.20) must have period 27t/q; in z, it follows that 
the only wave numbers, cr, corresponding to eigenfunctions having the required 
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period in z are those for which o 2 = p 2 cr 2 for some integer p. For each 
p(p = 1,2, * * • ) the eigenvalue problem (4.21) has an infinite sequence of real, 
nontrivial solutions 

(fb<£i,<£ 2 ) = (M P q,0i Pq ^D> P = 1,2, • • • ; q = ±1,±2, • • • . 

Since (— n ,— is a solution of (4.21) whenever is a solution of 

(4.21), we may order the indices so that 

(4.23) </>f ( ~ q) = -</>f q , </>| (_q) = <f>p, // p( _ q) = -/x pq , and 0 < n P i < n p2 < 
Using this notation, we see from our assumption (4.22) on a that 

(4.24) Ho = min ^ p i = M Po i 

P 

so that ^ pq > Ho (p,q) ± (Po>!)> Q > 1- 

The above discussion of the underlying problem (4.21) shows that the full 
eigenvalue problem (4.19) in H has the solutions 

(4.25) \ = H pq and u = #* q , P = 1,2, * * * ; q = ±1,±2, • * ■ , 
where 

(4.26) i[P q = (^f q cospaz, <j>P cospaz, ^fVmpaz) 

with <f>p determined by (4.20c). It follows as in [12, Appendix] that the eigen- 
functions {^P q } may be assumed orthonormal in H, after rescaling by con- 
stants depending upon p and q, i.e., 

(r>,<n - v<»- 


(4.27) 
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where <5 ik is the usual Kronecker delta symbol. 

The following lemma summarizes some of the basic facts for the linearized 
problem (4.19). The compactness properties are essentially known (e.g., see [10; 
11]) while the characterization (4.28) follows easily from (4.17). 

Lemma 4.1. (i) The linear operator L: H —*■ H defined in (4.16) is 
self-adjoint and compact and its characteristic values and eigenfunctions are 
given by (4.25). The eigenfunctions {t^ q } satisfy (4.27) and are complete in H. 

(ii) The linear operator M: H — ► H is compact and its adjoint, M*, is 
characterized by 

(4.28) (M*v,w) = 2 J yv^, v,w € H . 

0 

We show next that hypotheses (Hi) through (H4) of Section 1 are satisfied 
so that we can make use of the structure parameter approach developed in Sec- 
tion 2. To minimize the calculations, we shall make repeated use of the results 
obtained in [12; 13]. This may be done in most cases simply by replacing <f> 3 
and (f > 4 in [12] by (f> 2 and <f> 1 of the present paper, respectively. 

Since defined in (4.24) is a simple eigenvalue of (4.21) and /u 0 < 
for (p,q) # (p 0 ,l), q > 1, H o is simple and also the smallest positive charac- 
teristic value of L in H. The associated null space, N, of I— /%L is spanned 
by = ^ ol and TV"*" is spanned by (p,q) ^ (Po4)- Thus, making use 

of part (i) of Lemma 4.1 above, we see that (Hi) is satisfied. If M is defined as 
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in (4.17), then (H2) can be verified by using part (ii) of Lemma 4.1 above and 
part (ii) of Lemma 3.1 in [12]. The form of X c = X c (7) in (H3) may be derived 
as in Lemma 3.2 of [12]. In fact, the characteristic values of L 7 = L— 7M are 
determined by the problem obtained from (4.21) by replacing (4.21a) by 

(4.29) DV 2 — /zcr^l— : 2'-yy)0i = 0. 

Thus, X c is simple and is also real as an eigenvalue of the problem (4.21) with 
(4.21a) replaced by (4.29) and <7 set equal to a 0 (e.g., see [8]). Finally, setting 
G(w) = 0 in (H4) and recalling the definition of $ in (4.18), one sees that F in 
(4.15) is generated by the bounded bilinear operator (see also [10; 11]). 
Since it follows as in part (iv) of Lemma 3.1 in [12] that H X N — ► n\ we 
see that (H4) is also satisfied. 

To determine the coefficients a and c in the bifurcation equation (2.4), 
we note first of all that $ satisfies, in addition, the condition (1.3) in Remark 

1.1 (see part (iii) of Lemma 3.1 in [12]). Thus, we may make use of parts (v) 
and (vi) in Lemma 3.1 of [12] to show that a = 0 always and c > 0 in essen- 
tially all cases. Moreover, by making use of the invariance of equations (4.9) 
under the translation z — ► z + 7r /a, one can show that, for each r and 7, the 
remainder term r in (2.5) is odd in 0 (see the last part of Appendix B in [19] 
and also the proof of Lemma 3.2 in [13]). Thus, r has the special form given in 
(3.1), where, in addition, s is even in 0 . 

In view of the above discussion one can now use Theorem 3.1 and Corollary 

3.1 with 0 < 7 < 7 0 to determine nontrivial solutions of the operator equation 
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( f ) and, hence, roll-like solutions of the disturbance equations (4.9) satisfying 
the boundary conditions (4.10). Thus, for each fixed 7 satisfying 0 < 7 < 7 0 , 
where 7 0 is sufficiently small, rotating Couette-Poiseuille flow is stable up to 
X c (7) at which point it loses stability to a supercritical, stable roll-type solution 
given by (2.7) with /?*=/?, r* given by (3.4) with a = 0, and Vb = given 
by (4.26). 

As we now show, the above result implies that the superposition of a 
Poiseuille flow on a rotating Couette channel flow is, in general, destabilizing. 
Recall that X c (7) is given by an expression such as (1.2) with 
b = (MKMt/b,*/fi)* If b > 0, then, for 7 sufficiently small, X c (7) < X c (0) = /i 0 , 
where f.i 0 is the critical eigenvalue of the linearized problem for rotating 
Couette flow; to see that b > 0 here, one can make use of the results in [3; 
§71(d)j (see also [6, p. 98]) for narrow-gap Taylor problems with the parameter 
(l— yu)/(l+/i) in [3] replaced by 7 in (4.8). Thus, for a given value of the swirl 


S, 0 < S < — , the addition of a sufficiently small component of Poiseuille flow 

2 


to a basic rotating Couette channel flow always leads to roll-type solutions at 
values of X that are greater than X c (7) but less than the critical eigenvalue, 
/i 0 , at which rotating Couette flow loses stability. This type of result was con- 
jectured in [17] for non-rotating combined Couette-Poiseuille flow on the basis 
of numerical calculations for the linearized Navier-Stokes equation. 
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